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I .  INTRODUCTION 


The  interference  by  molecular  contaminants  in  the  environment  to 
the  proper  performance  of  space-based  infrared  telescopes  mounted  on  the 
Shuttle  Orbiter  has  been  discussed  by  a  number  of  authors.^  Contaminants 
include  natural  atmospheric  species,  primarily  oxygen  (0)  atoms  at  Shuttle 
altitudes,  as  well  as  molecules  released  by  the  Shuttle.  These  include  water 
(H2O),  nitrogen  (N2),  oxygen  (O2 ) *  carbon  monoxide  (CO),  carbon  dioxide  (CO2), 
other  trace  molecular  species,  and  hydrogen  (H2)  and  helium  (He),  although 
these  latter  two  will  not  normally  interfere  with  the  telescope's  operation. 
Potential  interference  effects  include  background  radiation  from  excited 
molecules  (and  possibly  from  interactions  of  0  atoms  with  surfaces,  i.e., 
"Shuttle  glow”),  as  well  as  degradation  in  optical  performance  due  to  con¬ 
densation  on  sensitive  optical  components,  specifically  the  cryogenically 
cooled  primary  mirror.  The  objective  of  the  work  reported  here  was  to  analyze 
the  degree  to  which  a  purge  flow  of  incondensible  He  gas  could  protect  against 
such  contamination  of  the  primary  mirror  of  the  AFGL  CIRRIS  (Air  Force 
Geophysics  Laboratory  Cryogenic  Infrared  Radiance  Instrument  for  Shuttle) 
telescope . 

Either  He  or  neon  (Ne)  could  potentially  be  used  as  the  incon¬ 
densible  purge  gas.  Although  Ne  should  have  better  efficiency  by  virtue  of 
its  higher  mass,  He  has  been  the  gas  chosen  in  virtually  all  discussions, 
analyses,  and  implementations  of  the  purge  concept  to  date.  One  practical 
reason  in  the  case  of  CIRRIS  is  the  use  of  a  liquid  He  cryogen  so  that  the  He 
gas  is  already  available  at  flow  rates3  m  up  to  0.02  g/sec.  Early  analyses  of 
the  He  purge  concept,  reviewed  in  Refs.  2  through  4,  demonstrated  its  feasi¬ 
bility,  but  were  too  approximate  to  calculate  or  measure  its  efficiency  reli¬ 
ably.  These  early  analyses  had  to  assume  a  simple  hard-sphere  differential 


aActually,  the  rate  of  gaseous  He  boil-off  from  the  cryogen  is  considerably 
faster,  but  only  10%  of  that  gas  is  used  in  the  He  purge  flow  in  the  current 
CIRRIS  design. 
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scattering  cross  section.  Accurate  differential  cross  sections  o( 6)  for  the 
scattering  of  He  atoms  by  0  atoms  reported  in  Ref.  3  show  properties  very 
different  from  the  simple  hard-sphere  approximations.  References  A  and  5 
provide  the  only  analyses  that  use  the  accurate  He-0  scattering  cross  sections 
calculated  in  Ref.  3.  Because  it  formulates  the  problem  in  terms  of  scattered 
0  atoms  that  miss  the  mirror  rather  than  those  that  hit  it,  Ref.  A  avoids  the 
unrealistic  limit  of  only  partial  purge  efficiency,  even  in  the  limit  of  infi¬ 
nite  He  density  which  appears  in  Ref.  5.  However,  both  calculations  employ  a 
single-collision  approximation  and  are,  consequently,  invalid  at  high  purge- 
flow  conditions. 

The  goal  of  the  present  study  is  to  calculate  the  He  purge-flow 
efficiency  without  the  limitations  of  the  previous  studies.  Specifically,  a 
Monte  Carlo  calculation  is  made  of  the  efficiency  of  the  He  purge,  as  a  func- 
tion  of  He  gas  flow  rate  m,  in  protecting  the  primary  mirror.  The  calculation 
is  valid  for  any  arbitrary  number  of  collisions  between  He  atoms  and  an 
incoming  contaminant  atom  or  molecule.  The  best  available  information  is  used 
to  calculate  highly  reliable  scattering  cross-section  functions,  and  purge 
protection  is  calculated  for  impinging  contaminant  molecules  released  from  the 
Shuttle  Orbiter,  as  well  as  for  impinging  atmospheric  0  atoms. 


II.  PROCEDURE 


Calculations  were  done  for  the  CIRRIS  I  and  CIRRIS  IA  tele¬ 
scopes.  The  geometry,  taken  from  Ref.  4,  is  shown  in  Fig.  1  for  CIRRIS  I. 

The  computer  code  developed  in  this  study  can  treat  the  case  of  Re  injection 
within  the  telescope,  or  He  ejection  from  orifices  in  front  of  the  telescope 
that  make  any  angle  u>  with  the  sensor  axis.  All  calculations  reported  here 
were  done  for  the  purge  apparatus  used  on  CIRRIS  I  and  CIRRIS  IA,  an  injector 
ring,  with  u>  =  30°,  about  the  circumference  of  the  sensor  entrance  aperture. 

A  plane  wave  of  incident  contaminant  molecules  of  mass  mc  approaches  the 
sensor  with  speed  vc  at  an  angle  y  to  the  sensor  axis.  The  purpose  of  the  He 
purge  gas  is  to  deflect  incoming  contaminants  that  would  otherwise  hit  the 
mirror  at  the  rear  of  the  sensor.  Detailed  analyses  are  presented  for  and 
H2O,  the  two  most  troublesome  contaminant  species.  Contaminant  molecules  may 
collide  with  He  gas  within  the  telescope  and  in  the  "plume,"  i.e.,  in  front  of 
the  telescope.  However,  He  density  falls  off  rapidly  in  the  plume  with 
increasing  distance  from  the  sensor.  In  these  calculations,  collisions  were 
considered  within  a  hemispherical  plume  region  of  100  cm  radius  centered  on 
the  sensor  aperture. 

The  temperature  of  the  CIRRIS  primary  mirror  is  20  K;  an  analvsis 
presented  later  in  this  report  indicates  that  even  will  cond<>nse  at  rM  = 
temperature.  The  temperature  of  the  walls  of  che  CIRRIS  telescope  varies  from 
about  20  K  near  the  primary  mirror  to  about  40  K  near  the  sensor  aperture. 
Contaminants  such  as  N2 ,  which  condense  only  at  very  low  temperature,  will  not 
condense  along  the  telescope  wall  except  at  the  back  near  the  primary  mirror. 
On  the  other  hand,  the  vapor  pressure  of  HjO  is  negligible  even  at  40  K. 
Accordingly,  it  was  assumed  that  an  impinging  H2O  molecule  sticks  to  the  first 
surface  it  hits,  mirror  or  wall.  This  same  assumption  has  been  made  in  all 
previous  treatments  of  the  He  purge  concept.  It  should  be  noted  that  a  record 
of  the  profile  of  contaminant  deposition  along  the  sensor  wall  was  kept  for 
all  calculations  reported  here.  Consequently,  it  would  be  possible  in  the 


future  to  examine  sensitivity  to  this  assumption  of  a  sticking  coefficient  oc 
unity  by  allowing  contaminants  to  reflect,  after  partial  accommodation,  from 
various  wall  positions. 

A.  HELIUM  FLOW  FIELD 

The  spatial  dependence  of  the  He  number  density  and  velocity  dis¬ 
tributions  are  necessary  input  functions  for  the  calculations  of  the  He  purge 
efficiency.  These  functions  were  modeled  by  assuming  thermal  effusion  of  He 
from  the  ejection  ring  shown  in  Fig.  1.  Perturbations  of  the  flow  field  that 
resulted  from  He-He  collisions  outside  of  the  ejection  ring  were  ignored. 
Possible  perturbations  of  the  He  flow  field,  which  resulted  from  collisions 
with  contaminants  or  atmospheric  species,  were  also  ignored.  This  latter  con¬ 
dition  should  be  an  excellent  approximation  at  the  Shuttle  Orbiter  altitude, 
380  km,  planned  by  CIRRIS,  although  it  could  break  down  at  much  lower  orbital 
altitudes . 

The  temperature  of  the  He  in  the  ejection  ring  in  CIRRIS  IA  may 
vary  from  ~  220  to  273  K.  This  relatively  high  temperature  resulcs  because 
boil-off  from  the  liquid-He  cryogen  is  used  as  a  coolant  for  an  electronics 
box  before  it  enters  the  ejector  ring.  Most  of  the  calculations  were  carried 
out  for  a  He  ejection  temperature  T^g  of  273  K.  However,  some  calculations 
were  carried  out  for  a  much  lower  ejection  temperature,  T^e  =  20  K,  In  order 
to  examine  the  sensitivity  of  He  purge  efficiency  to  ejection  temperature. 

At  any  point  in  the  plume,  the  He  number  density  was  computed  by 
summing  the  effusion  contribution  from  sources  on  the  circumference  of  the 
ejection  ring  and  effusion  from  the  sensor  entrance  aperture.  The  He  number 
density  within  the  sensor,  assumed  to  be  a  constant  within  the  sensor  volume, 
was  calculated  by  equating  rates  of  effusion  into  and  out  of  the  sensor  aper¬ 
ture.  In  order  to  save  computer  time  and  render  the  purge  problem  computa¬ 
tionally  tractable,  all  effusing  He  was  assumed  to  move  at  an  average  thermal 
speed  of 

vHe  -  (8  kT/Tnn^)172  (1) 


Within  the  plume,  the  distribution  in  directions  of  He  motion  was  obtained  by 
summing  all  effusion  contributions.  As  noted  above,  there  is  a  distribution 
in  temperatures  within  the  sensor.  In  treating  the  He  flow  field,  however, 
this  was  treated  as  one  temperature,  Ty.  The  distribution  in  speeds  within 
the  plume  could  be  bimodal,  corresponding  to  the  sensor  and  ejection-ring 
temperatures.  The  speed  distribution  anywhere  within  the  sensor  consisted  of 
He  moving  at  v^g  from  Eq.  (1),  with  a  random  distribution  in  directions. 

Figure  2  shows  a  comparison  between  the  predictions  of  this  number 
density  model  and  measurements  reported  in  Ref.  4  for  a  0.0405-g/sec  ring 
ejection  of  300  K  argon  (Ar)  gas  in  a  1/1.78-scale  CIRRIS  I  model.  The  agree¬ 
ment  at  large  distances  is  encouraging  because  model  and  measurements  are 
steady-state  results  for  the  same  mass  flow  rate  m,  so  that  the  integrals  of 
the  net  outward  flux,  over  any  bound  surface  containing  the  aperture,  have  to 
be  equal.  The  discrepancy  at  smaller  distances  might  be  due  to  noneffusive 
real  source  flow,  purge-purge  collisions  in  the  real  flow,  as  well  as  experi¬ 
mental  error.  The  approximation  of  effusive  flow  from  the  orifices  in  the 
ejector  ring  would  result  in  the  largest  uncertainty  in  He  densities  in  this 
region,  i.e.,  within  and  immediately  in  front  of  the  sensor  aperture.  In 
order  to  examine  the  sensitivity  of  purge  efficiency  results  to  He  flow-field 
uncertainties,  calculations  were  actually  carried  out  as  a  function  of  the 
steady-state  density,  denoted  N^g ,  of  He  within  the  sensor.  Additional 
sensitivity  studies  could  be  carried  out  by  examining  noneffusive  source  flow 
velocity  distributions.  However,  this  was  not  pursued  because  the  effect 
should  be  very  small  in  comparison  to  uncertainties  introduced  by  the  approx¬ 
imate  treatment  of  the  CIRRIS  IA  geometry,  which  is  discussed  later. 

B.  COLLISION  CROSS  SECTIONS 

A  contaminant  or  atmospheric  species  changes  speed,  usually 
slowing  down,  each  time  it  collides  with  a  He  atom.  The  collision  cross 
sections  are  needed  for  a  wide  range  of  relative  collision  speeds, 
v  »  | v  -  v  |,  in  order  to  follow  many  collisions  between  incoming 
contaminants  and  He  purge  atoms.  Differential  and  total  cross  sections 
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Measured  (Ref.  4)  and  Modeled  (see  text)  Centerline  Number 
Densities  for  Subscale  Ar  Ring  Injection.  Z  >  0  measures 
distance  in  front  of  sensor  aperture.  Results  for  Z  <  0 
indicate  values  within  the  sensor. 
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a(v)  =  2 tt  /  a(9,v)  sin6d0  (2) 
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for  the  scattering  of  0  atoms  from  He  atoms  were  calculated  in  Ref.  3  by 
solving  the  quantum  mechanical  Schrodinger  equation  for  the  best  information 
available  in  the  literature  on  the  central  force  interaction  potential  V(r). 
Sensitivity  analyses  reported  in  Ref.  3  illustrate  inconsequential  errors 
in  a(8,v)  and  o(v)  that  result  from  the  assumption  of  a  single  central  force 
potential  and  from  possible  errors  in  the  potential  function.  The  literature 
indicates  that  the  interaction  potentials  between  He  and  the  molecules  CO,^ 
N2,^  and  are  only  weakly  anisotropic.  To  an  excellent  approximation, 

they  can  be  analyzed  as  isotropic  systems;  i.e.,  the  interaction  potential  is 
assumed  to  depend  only  on  r  and  not  on  the  orientation  (central  field  approxi¬ 
mation)  of  the  molecule.  The  total  and  differential  scattering  cross  sections 
for  thes®  "isotropic  cases”  (i.e.,  He  +  Nj ,  He  +  CO,  and  He  +  ^0)  were  calcu¬ 
lated  quantum  mechanically  by  means  of  the  standard  techniques  described  in 
Ref.  3.  Results  for  the  total  scattering  cross  section  are  shown  in  Fig.  3. 
Appendix  A  provides  plots  of  the  differential  cross  sections  and  describes  the 
interaction  potentials  used. 

Q 

The  He  +  CO2  system  is  strongly  anisotropic,  and  the  "central 
field  approximation"  could  not  be  used  in  the  present  study.  The  differential 
cross  section  in  this  case  was  calculated  by  means  of  the  infinite  order 
sudden  (IOS)  approximation.^  The  interaction  potential  V(r,ct)  is  a  function 
of  the  separation  distance  r  and  the  angle  a  between  the  line  of  center  and 
the  axis  of  the  CC>2  molecule.  In  the  IOS  approximation,  a  is  regarded  as  a 
parameter  and  not  a  dynamical  variable.  For  a  fixed  value  of  the  parameter 
a,  the  potential  is  a  function  only  of  r  and  is  thus  "isotropic."  This 
"isotronic"  scattering  problem,  for  a  fixed  a,  is  then  solved  by  the  partial- 
wave  method  to  obtain  a  differential  cross  section  o(0,a),  which  is  also  a 
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function  of  the  parameter  a.  It  has  been  shown  elsewhere  that  the  desired 
"total"  (elastic  plus  inelastic)  differential  scattering  cross  section  is,  to 
an  excellent  approximation,  given  simply  by  the  orientation  average  over  a: 

1  * 

o(9)  ■  7  /  0  (0,o)  sin(a)  da  (3) 

0 
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In  the  present  work,  the  above  integral  was  calculated  numerically  by  means  of 
the  trapezoidal  rule,  with  an  integration  spacing  of  Act  *  22.5°.  By  making 
use  of  symmetry,  the  integration  range  was  reduced  from  0  to  it  to  0  to  tt/2  . 
Thus,  four  separate  partial-wave  calculations  to  compute  0(6,0^)  for  = 
22.5°,  45°,  67.5°,  and  90°  (a  *  0°  need  not  be  calculated  because  of  the 
sina  factor)  were  calculated  to  obtain  each  differential  scattering  cross- 
section  curve  for  He  +  CO2 ,  shown  in  Appendix  A.  Again,  the  total  scattering 
cross-section  results  are  shown  in  Fig.  3. 

C.  MONTE  CARLO  CALCULATIONS 

A  Monte  Carlo  calculation  was  made  for  contaminant  molecules 
impinging  at  angle  y  with  speed  vc  (see  Fig.  1).  The  transmission 

number  of  contaminants  striking  mirror  with  He  flow  . 

number  of  contaminants  entering  sensor  without  He  flow 


is  the  quantity  of  interest  in  analyzing  readily  condensible  contaminants, 
such  as  H2O.  For  contaminants  that  condense  only  at  the  primary  mirror,  viz 
N2,  the  quantity  of  interest  is  the  sensor  acceptance 

n  =  number  of  contaminants  entering  sensor  with  He  flow  .... 

1  number  of  contaminants  entering  sensor  without  He  flow 

The  Monte  Carlo  technique  is  a  well-known  procedure^®-^  wherein  unspecified 
parameters  are  chosen  at  random,  and  a  large  number  of  individual  trajectories 
of  contaminants  are  followed  as  they  transverse  the  He  flow  field,  in  order  to 
provide  a  statistically  reliable  result. 

Specifically,  the  calculation  placed  a  control  surface  normal  to 
the  impingement  direction  in  front  of  the  sensor  at  a  point  where  the  He  plume 
density  was  ~  IZ  of  its  centerline  value  at  the  sensor  aperture.  The  area  of 
the  control  surface  through  which  the  contaminants  impinged  was  typically 
twice  that  of  the  sensor  aperture,  in  order  to  allow  for  scattering  of  con¬ 
taminants  into,  as  well  as  out  of,  the  sensor  aperture.  A  record  was  kept 
for  each  calculation  of  the  number  of  trajectories  where  the  He  actually 


scattered  an  impinging  contaminant  into  the  sensor  aperture.  This  was  always 
a  small  fraction  of  the  trajectories,  a  fact  confirming  that  the  control 
surface  area  was  sufficient  to  account  for  this  effect.  This  was  further 
verified  in  a  few  calculations  where  other  control  areas  were  used. 

The  spatial  position  of  the  contaminant,  as  it  crossed  normal  to 
the  control  surface,  was  chosen  at  random.  The  contaminant  was  allowed  to  move 
along  a  straight-line  trajectory  s  for  a  distance  £  before  a  collision  took 
place.  The  free  path  length  i  was  calculated  from  the  randomly  chosen 
variable,  0  <  x  <  1,  as 

1  1 

—  /  <vo(v)>  n(r)ds  *  -Inx  (6) 


<vo(v)>  =  /  vo(v)p(v)dv  (7) 

0 

is  the  numerical  average  of  the  total  collision  rate  constant  over  the 
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probability  density  distribution  in  relative  collision  speed,  viz,  J0  p  (v) 
dv  ”  1 .  The  quantity  p(v)  has  two  forms  —  one  for  outside  the  sensor,  where 
the  He  velocity  field  is  determined  by  effusive  flow  from  the  ejector  ring  and 
sensor  aperture,  and  one  for  inside  the  sensor,  where  the  He  velocity  field  is 
random  in  direction.  In  both  cases,  p(v)  is  for  the  mean-speed  distribution 
approximation  f(v^e)  ”  <S(v^e  “  vHe^’  the  resu^tant  straight-line  trajec¬ 

tory  corresponded  to  hitting  the  sensor  wall,  hitting  the  mirror,  or  exiting 
the  He  plume  to  free  space,  the  result  was  recorded  and  the  motion  of  the  next 
incoming  contaminant  was  followed. 

If  the  solution  for  £  resulted  in  a  position  in  the  He  plume  or 
In  the  He  gas  within  the  sensor,  a  collision  had  taken  place.  After  the 
collision,  the  contaminant  left  this  position  on  a  new  straight-line  trajec¬ 
tory  with  a  new  velocity  v^,  i.e.,  in  a  new  direction  and  at  a  new  speed,  and 
the  procedure  was  repeated  to  find  a  new  trajectory  length.  First,  however, 

the  new  contaminant  velocity  v'  was  calculated  from  the  old  known  contaminant 
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velocity  vc  by  invoking  conservation  of  energy  and  linear  momentum  in  the 
collision.  For  this  calculation  it  was  necessary  to  specify  the  precollision 
direction  of  motion  of  the  He  atom  involved  in  the  collision,  as  well  as  the 
center-of-mass  scattering  angles,  0  and  produced  by  the  collision.  A  value 
for  the  azimuthal  angle  $  was  chosen  randomly,  i.e., 

^  =  2mc  (8) 

The  polar  scattering  angle  9  was  chosen  by  randomly  selecting  a  cumulative 

scattering  probability,  i.e.,  by  solving 

2  0 

P  (9)  -  —  /  o(  0'  )sin0'd  9'  -  x  (9) 

c  <3  JQ 

for  0  when  x,  as  usual,  is  randomly  chosen  between  0  and  1.  Rigorously,  the 
He  direction  should  be  chosen  randomly.  This  was  done  for  collisions  within 
the  sensor.  For  collisions  within  the  plume,  however,  the  He  direction  was 
taken  as  the  most  probable  one  for  that  particular  spatial  position  in  order 
to  reduce  computation  time.  The  effect  of  this  approximation  on  the  compu¬ 
tation  should  be  negligible. 

In  this  manner,  trajectories  of  a  large  number  of  impinging  con¬ 
taminants  were  followed  to  arrive  at  statistically  significant  results  for  T 
and  n.  Each  trajectory  was  made  up  of  an  arbitrary  number  of  straight-line 
trajectories,  with  intervals  between  collisions  determined  randomly.  After 
each  collision,  the  new  contaminant  velocity  was  determined  randomly,  consis¬ 
tent  with  conservation  of  linear  momentum  and  energy. 

Equation  (9)  indicates  that  the  differential  cross  section  only 
enters  the  Monte  Carlo  calculation  of  r  through  the  cumulative  scattering 
probability  distribution  Pc(0).  Figure  4  shows  plots  of  this  function  for  He 
+  CO2  for  three  different  collision  velocities.  For  all  of  the  Monte  Carlo 
calculations  reported  here,  Pc( 0)  was  fit  to 
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Fig.  4.  Plots  of  Pc  -  (2* /o)  /0  o(d’)  sin0'd0’  for  He  +  C02  for  Three 
Relative  Collision  Speeds  v 
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As  n  goes  to  unity,  the  scattering  would  approach  isotropic,  hard-sphere 
scattering.  As  n  goes  to  infinity,  the  scattering  is  more  and  more  closely 
confined  to  a  region  around  0°.  These  two  limiting  scattering  behaviors, 
isotropic  and  pure  forward,  are  also  shown  in  Fig.  4.  Figure  5  shows  the 
dependence  of  n  on  collision  speed  for  He  +  HjO,  He  +  N2 ,  and  He  +  CC>2 , 
inferred  by  fitting  the  data  in  the  Appendix  to  Eq.  (10).  The  fact  that  n  is 
much  greater  than  unity  over  the  range  of  collision  speeds  of  interest  here 
illustrates  the  failure  of  the  hard-sphere  scattering  approximation  for  this 
particular  problem. 
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III.  ESTIMATED  CONTAMINATION  LEVEL 


The  characteristics  of  the  expected  contaminants  are  necessary 
input  conditions  for  the  He-purge  transmission  calculations.  At  an  orbital 
altitude  of  380  km  and  a  velocity  of  vq  *  7.7  x  10^  cm/sec,  0  atoms  will  be 
the  most  abundant  atmospheric  species.  For  standard  atmospheric  conditions, 

Q 

they  will  be  present  at  a  density  of  n  =  1.7  x  10  cm  ,  corresponding  to  a 
flux  of  F  »  Vgti  =  1.4  x  10 14  cm-^  sec-^  .  Rcwever  *  they  could  reach  the  CIRRIS 
primary  mirror  only  for  a  telescope  line  of  sight  that  makes  a  small  angle 
0  with  the  direction  of  the  orbital  velocity  vQ,  whereas  the  CIRRIS  IA  mission 
plans  specify^  that  0  >  90°.  For  these  large  values  of  0,  the  0  atoms  inci¬ 
dent  along  the  telescope  axis  will  be  negligible.  Indeed,  it  can  be  shown^ 
that  the  incident  flux  decreases  very  rapidly  with  increasing  0  as 

f  *  K  {a(2a2  +  3)  [erf (a)  +  1]  exp(-c2sin2 0)  + 


(2/tt1/2)  (a2  +  1)  exp(-c2)  } 


1  2  —2  —1  —1 

where  K  -  1.6  *  10  4  0  atoms  cm  ^  sr  sec  ,  a  =  c  cos  0,  and  c  *  7.0  is  the 
ratio  of  vQ  to  the  most  probable  0  atom  thermal  speed.  At  M  0°,  f  ■  2.2  x 
10  0  atoms  cm  sr  sec  ,  whereas  f  <  5  x  10  0  atoms  era  sr  sec  for 

any  0  between  90°  and  180°,  i.e.,  attenuated  by  more  than  23  orders  of 
magnitude.  It  should  also  be  noted  that,  by  virtue  of  the  comparable  masses 
and  similar  cross-section  behaviors  indicated  in  the  Appendix,  He  purge 
results  that  are  calculated  here  for  1^0  would  be  approximately  valid  for  0 
atoms  as  well. 

Reference  16  reports  results  from  a  mass  spectrometer  flown  on 
the  Shuttle  Orbiter.  The  inlet  of  the  mass  spectrometer  was  collimated  to 
a  narrow  field  of  view,  0.1  sr,  of  the  space  environment  so  that  molecules 
outgassing  from  Orbiter  surfaces  could  not  reach  it  by  a  straight-line  trajec¬ 
tory.  Contaminants  entering  the  CIRRIS  telescope  are  under  the  same  geometric 
constraint;  contaminants  enter  either  instrument  only  by  virtue  of  some  other 
transport  mechanism,  such  as  self-scattering  or  scattering  from  atmospheric 
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0  atoms.  Detected  contaminants  reported  in  Ref.  16  include  H? ,  He,  diborane, 
(B9H^),  methane  (CH^),  N-> ,  CO,  0,  monomethyl  hydrazine  (CH^NSH-j),  H  ?0,  C0? , 
nitric  oxide  (NO),  0^,  hydrogen  chloride  (HC1),  and  propane  (C3H^).  With 
the  possible  exception  of  N£,  which  is  discussed  later,  H2O  was  the  dominant 
species.  This  same  predominance  of  contaminant  H^O  was  also  observed^  by  a 
different  mass  spectrometer  instrument  on  one  Shuttle  flight;  in  this  case, 
however,  outgassed  molecules  could  reach  the  mass  spectrometer  directly  by 
collisionless,  straight-line  trajectories. 

Table  1  lists  the  1^0  contamination  levels  reported  in  Ref.  16  for 
Space  Transportation  System  (STS)  flights  3  and  4.  The  source  of  this  H2O  was 
outgassing  from  Shuttle  Orbiter  surfaces.  Typically,  outgassing  rates  are 
high  immediately  after  launch  and  decay  over  a  period  of  many  hours  to  a 
steady  value.  Values  listed  in  Table  1  are  approximate  100-hr  averages  for 
the  sensor  perpendicular  to  the  wind  direction,  i.e.,  for  g  =  9D°.  There 
have  also  been  a  number  of  predictions  of  the  expected  Orbiter  contaminant 
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environment.  Two  representative  predicted  values  ’  are  also  included  in 
Table  1  to  indicate  that  the  measured  values  conformed  to  reasonable 
expectations. 

Thus,  CIRRIS  IA  will  probably  have  to  contend  with  a  contaminant 
HtO  environment  similar  to  that  present  on  STS  flight  3,  i.e.,  ~  3  x  10^  H,0 

_  O  _  1  _  1 

molecules  cm  sr  1  sec  .  However,  there  is  no  doubt  that  the  contaminant 
H9O  environment  was  worse  on  STS-4  because  of  the  hail  and  rain  experienced  by 
Columbia  during  the  STS-4  prelaunch  preparation.  It  is  therefore  prudent  to 

19  _  9  _  1  _  1 

take  the  STS-4  environment,  ~  7  x  10l  H2O  molecules  cm  “  sr  sec  ,  as  the 

postulated  worst-case  environment  for  the  CIRRIS  IA  mission. 

Transmission  through  the  He  purge  flow  depends  on  the  speed  vc  of 
the  incident  contaminant.  This,  in  turn,  depends  on  the  mechanism  that  trans¬ 
ports  the  contaminant  from  its  source  to  the  sensor  aperture.  For  a  contami¬ 
nant  of  mass  mc  that  escapes  from  the  Shuttle  Orbiter  and  scatters  off  an 
atmospheric  0  atom  of  mass  mg  back  toward  the  sensor 
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HpO  Flux  (cm-  sr-  sec-') 


Outgasslng  -  STS-3 


Ref. 

16 

Mass  Spectrometer  Measurement 
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Ref. 

18 

Calculation 
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Ref. 

19 

Calculation 

2 

X  1010 

Ref.  16 


Ref.  17 
Ref.  18 
Ref.  19 


Outgasslng  -  STS-4 
Mass  Spectrometer  Measurement 

Vernier  Engine  Firing 

Mass  Spectrometer  Measurement* 

Calculation 

Calculation 


7  x  10J 


0.4  -  4  x  10 

1  x  1014 
1  x  1013 


aThese  are  approximate  values  for  a  sensor  line  of  sight  perpendicular 
to  the  wind.  Actual  values  reported  in  Refs.  16  through  19  were  scaled 
to  those  at  380  km  by  assuming  contaminant  return  flux  to  be  proportional 
to  atmospheric  density. 

^Calculated  by  assuming  a  2n-steradian  field  of  acceptance  for  the  mass 
spectrometer . 
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for  0°  <  3  <  90°.  This  expression  fails,  of  course,  near  B  =  90°,  where  it 
predicts  that  vc  goes  to  zero.  The  contaminant  flux  from  this  mechanism  would 
also  be  3  dependent,  peaking  at  6=0°  and  dropping  to  zero  at  3  =  90°. 
Reference  16  does  provide  some  qualitative  data  on  this  3  dependence.  The 
values  in  Table  1  from  Ref.  16  refer  to  some  3  s  90°.  Data  in  Ref.  16 
actually  indicate  that  the  H20  return  flux  varies  by  a  factor  of  ~  80  between 
3=0°  and  90°.  Calculations  reported  in  Ref.  19  are  in  reasonable  agreement, 
giving  a  factor  of  ~  25  for  this  variation  in  return  flux.  The  variation 

in  H2O  return  flux  with  3  is  probably  a  sharp  function  of  3-  In  the  absence 
of  any  better  information,  and  in  order  to  make  a  conservative  estimate  of  the 
H2O  return  flux,  however,  the  H2O  return  flux  used  to  estimate  1^0  film  build¬ 
up  on  the  CIRRIS  IA  primary  mirror  is  taken  to  be 

fH  0<8)  =  f h  0(S  =  90<>)  [1  +  (?9  cosB)l’  0°  <  B  <  90° 

and 

fH20(B)  =  fH20(B  =  9°0) >  8  >  90°  (13) 

where  (90°)  =  7  x  1012  ^0  molecules  cm*2  sr-1  sec-1  (STS-4  result). 

For  3  >  90°,  the  return  flux  is  due  to  self-scatter  of  two  or  more  collisions 
with  the  wind,  and  it  should  return  with  speeds  comparable  to  or  less  than 
thermal  speeds  of  outgassed  molecules.  Thus,  for  angles  3  >  90°,  vc  is  taken 
as  ~  6  x  10^  cm/sec ,  the  average  thermal  speed  of  at  300  K. 

The  CIRRIS  IA  mission  plan  currently  specifies  that  no  engine  in 
the  Orbiter  may  fire  while  the  telescope  is  uncapped.  This  includes  the  small 
25-lb-thrust  vernier  control  system  (VCS)  engines.  In  order  to  estimate  the 
consequences  if  this  requirement  were  relaxed,  the  contaminant  environment 
produced  by  these  engines  has  also  been  considered.  The  largest  source  of 
uncertainty  is  the  return  speed  of  molecules  that  are  produced  when  these 
engines  fire,  as  there  are  neither  measurements  nor  calculations  of  this 
parameter.  Along  the  axis  of  the  plume,  the  flow  velocity  is  ~  3.5  x  108 
cm/sec.  However,  molecules  that  return  to  the  Orbiter  bay  probably  originate 
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from  the  wings  of  the  plume  flow  field.  Here,  plume  flow  velocity  may  he  less 
than  the  flow  velocity  along  the  axis.  Again,  however,  there  are  neither  mea¬ 
surements  nor  calculations  of  velocities  in  this  region  of  the  plume  of  this 
engine.  In  any  case,  plume  molecules  can  return  to  the  Orbiter  Bay  because  of 
collisions  with  the  atmosphere  at  a  wide  range  of  speeds  vc  that  correspond  to 
the  range  in  possible  geometries  of  motor  location  and  thrust  direction  in 
relation  to  the  Orbiter's  attitude  and  the  sensor's  line  of  sight. 

For  H2O,  a  reasonable  range  of  return  speeds  is  v^  “  0.6  -  5  x  10^ 
cm/sec.  Reference  17  measured  density  spikes  of  ^0  and  N£ ,  but  not  of  CO  or 
CO2,  associated  with  firings  of  these  VCS  engines.  The  corresponding  return 
fluxes  are  listed  in  Table  1;  the  range  in  values  quoted  corresponds  to  the 
range  of  uncertainties  in  vc*  Recently,  Ehlers  reported  calculations  of  the 
return  flux  as  a  function  of  VCS  engine  location.  His  calculations  exhibit  a 
large  variation  in  return  flux,  by  approximately  a  factor  of  AO,  as  the  angle 
between  the  Orbiter's  velocity  and  the  instrument's  line  of  sight  varies 
between  10°  and  90°,  in  approximate  agreement  with  the  factor  estimated  in 
Table  1.  When  converted  to  an  altitude  of  380  km,  however,  the  absolute 
magnitude  of  his  calculated  return  fluxes  are  lower  by  approximately  a  factor 
of  ten  than  the  estimate  given  in  Table  1.  Here  again,  the  estimate  in 
Table  1  is  used  to  calculate  the  rate  of  cryodeposit  buildup  in  order  to 
provide  the  more  conservative,  worst-case  analysis.  However,  it  is  important 
to  note  that  the  observations  reported  in  Ref.  17  included  contributions  from 
the  vernier  engines  on  the  rear  of  the  Orbiter.  The  plumes  of  these  motors 
impact,  at  least  partially,  on  the  wings  of  the  Shuttle.  Thus,  it  is  quite 
possible  that  data  from  Ref.  17  represent  an  unrealistic  worst  case  for  the 
CIRRIS  mission  because  direct  scattering  of  plume  molecules  from  the  Orbiter 
wing  could  have  contributed  significantly  to  the  observations  reported  in  Ref. 
17,  whereas  it  could  not  contribute  return  flux  along  the  CIRRIS  line  of 
sight.  For  this  reason,  cryodeposit  buildup  rates  based  on  Ehlers'^®  return 
fluxes  are  also  discussed  later. 

In  addition  to  ^0,  N2  is  probably  also  an  important  contaminant 
molecule  in  the  vicinity  of  the  Orbiter.  However,  the  information  on  the 
return  flux  of  this  contaminant  is  very  limited.  Owing  to  an  instrumental 


background  at  mass  28,  Ref.  16  reported  no  data  on  N2  return  flux  except  at 
small  8  angles,  where  impingement  of  atmospheric  N2  is  important.  Reference 
17  reports  N7  and  H2O  densities  that  are  in  a  ratio  of  ~  0.02  in  the  absence 
of  a  VCS  firing,  and  ~  1  in  the  transient  spike  in  the  presence  of  a  VCS 
firing.  In  the  absence  of  a  VCS  firing,  it  is  not  clear  what  fraction  of 

either  of  these  N2  or  1^0  signals  corresponds  to  a  return  flux.  In  a  recent 
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paper  recommending  environment  estimates  for  the  Orbiter  bay,  Scialdone 
recommends  an  on-orbit  N2/H2O  ratio  of  ~  AO,  in  clear  disagreement  with 
results  measured  in  Ref.  17.  Thus,  in  the  absence  of  data  on  the  N2  return 
flux,  the  value  adopted  here  is  the  calculated  flux  reported  in  Ref.  19  for 
6  -  90°.  This  was  adjusted  to  a  corresponding  flux  at  the  CIRRIS  altitude  of 
380  km  by  correcting  linearly  for  differing  atmospheric  densities  to  give 

f„  ( 8)  ■  1  x  10 11  molecules  cm  2  sr  ^  sec  ^ ,  8  >  90  0  (lAa) 
N2  2 

Here,  v£  is  taken  as  A. 6  x  10^*  cm/sec ,  the  average  thermal  velocity  of  N2  at 
300  K.  This  should  represent  a  conservative  overestimate  of  the  actual 
return  flux,  because  Ref.  19  calculated  it  for  an  N2  release  corresponding  to 
tie  maximum  possible  atmospheric  cabin  leakage  rate  that  NASA  specifications 
permit.  Again,  the  8  dependence  of  this  N2  flux  is  assumed  to  have  the  forms 

f  (8)  -  f  <B  -  90°)  {1  +  A9  cos  6)  ]  ,  0°  <  6  <  90°  (14b) 

2  2 

Once  more,  this  is  a  conservative  overestimate,  with  the  6  *  0°/6  *  90°  ratio 
exceeding  that  calculated  in  Ref.  19  by  a  factor  of  ~  2. 

A.  CRYODEPOSIT  BUILDUP  RATE 

For  a  contaminant  like  H20,  which  sticks  to  the  sensor  walls,  film 

buildup  on  the  primary  mirror  is  due  solely  to  molecules  that  strike  the 

mirror  first.  Consider  H20  impinging  through  a  point  on  a  hemisphere  of 

radius  r  that  is  centered  on  the  sensor  aperature.  The  solid  angle  subtended 

2 

by  the  sensor  aperture  is  A  cosy/r  ,  where  A  is  the  aperture  area.  Then,  the 
number  of  H20  molecules  that  pass  through  a  small  area  on  this  hemisphere  and 
would  enter  the  sensor  in  the  absence  of  any  He  flow  is  given  by 


fj 

r 


f  ( Y ,  (A  cosy/?  )r^  sinYdYd<}> 


For  Che  case  considered  here,  where  the  areas  of  Che  mirror  and  the  sensor 

aperture  are  the  same,  the  average  number  of  t^O  molecules  condensing  per 
2 

second  per  cm  of  mirror  becomes 

2  it  tt/2 

C  -  /  /  t(y)  f(Y,4>)  sin  Yd  Yd  $  (15; 

0  0 

As  discussed  in  Table  2,  the  factor  of  cosy  is  included  in  the  t(y)  values 

that  are  calculated . 

Clearly,  f(Y,$)  can  be  calculated  from  the  f(8)  expressed  in 
Eqs.  (13)  and  (14).  In  view  of  the  approximate  nature  of  this  assumed 
costs  functional  dependence,  however,  the  angular  flux  dependence  was  ignored 
in  evaluating  Eq .  (15)  to  obtain 


C  -  f  t(0)  B 


(15a) 


where 


b.Jl  f 

B  T(0)  i 


t(y)  sinYd y 


Nevertheless,  the  value  of  f  in  Eq .  (15a)  will  be  calculated  for  the  partic¬ 
ular  value  of  the  angle  8  between  v^  and  the  sensor  line  of  sight.  Then,  the 
rate  of  contaminant  film  buildup,  in  cm/ sec ,  is  given  by 

t  -  CM/6  x  1023  p  (17) 

where  M  is  the  molecular  weight  and  p  is  the  density. 

As  will  be  discussed  later,  N2  represents  the  opposite  case  of  a 
molecular  contaminant  having  a  temperature-dependent  vapor  pressure  that 
prevents  condensation  on  the  warmer  sensor  walls  but  that  allows  condensation 
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on  the  primary  ~irror.  In  this  case,  a  steady-state  N2  density  develops  in 
the  sensor  according  to 

n  =2  fn(0)  B'/v  (18) 

ss 

where 

,  n/2 

B’’^0)  /  n(Y)  sinydy  (19) 

v  is  the  average  N2  thermal  speed  of  1  *  10^  cm/sec  at  20  K,  and  a  factor  of  2 
enters  because  N2  enters  only  one  end  of  the  telescope  but  is  removed  at  both 
ends.  In  this  case,  the  cryodeposit  growth  rate,  in  cm/sec,  is  given  by 

t  -  0.25  n  v  M/6  x  1023  p  (20) 

ss 

t  «  0 

O  ft  OO 

where  n  *  4.8  x  10  cm  J  is  the  number  density  corresponding  to  the  vapor 
SS 

pressure  of  N2  at  20  K,  the  temperature  of  the  primary  mirror. 

B.  SIGNIFICANT  CRYODEPOSIT  THICKNESS 

Clearly,  some  small  contamination  buildup  on  the  primary  mirror 
can  be  tolerated.  However,  the  critical  thickness  associated  with  the  onset 
of  significant  degradation  in  sensor  optical  performance  is  not  well  defined. 

A  careful  definition  of  this  critical  thickness  would  itself  require  a  study 
comparable  in  effort  to  the  present  analysis  of  He  purge  efficiency.  This 
section  is  restricted  to  a  review  of  the  very  limited  published  information  on 
this  question  and  to  a  rough  estimate  of  the  thickness  of  cryodeposit  films 
which  might  begin  to  affect  the  optical  performance  of  the  CIRRIS  telescope. 

Both  absorption  and  scattering  by  the  contaminant  film  can 
adversely  affect  optical  performance.  Absorption  is  well  understood.  The 
thickness  of  an  ice  deposit,  for  example,  should  be  less  than  ~  0.1  un 
(1000  A)  in  order  to  prevent  significant  absorption  of  radiation  within  the 


H2O  absorption  bands. However,  the  scatter  properties  of  condensed  films 
are  only  poorly  characterized.  For  the  SIRE  (Satellite  Infrared  Experiment) 
sensor,  which  had  off-axis  rejection  requirements  similar  to  those  of  CIRRIS, 
it  was  estimated  that  cryodeposit  thicknesses  up  to  2  pm  would  produce  insig¬ 
nificant  scatter. ^  A  similar  estimate  of  1  in  for  the  SIRTF  (Shuttle  Infrared 
Telescope  Facility)  is  given  in  Ref.  1.  However,  the  basis  of  these  estimates 
is  not  documented;  they  probably  refer  to  a  homogeneous  smooth  deposit. 
Reference  24  specifies  a  bidirectional  reflectance  distribution  function 

(BRDF)  scatter  performance  for  the  primary  mirror  of  CIRRIS  of  lO-^  sr-^ 

_2 

at  0  *  1°,  with  a  0  rolloff.  These  authors  mention,  but  do  not  reference, 
measurements  that  indicate  that  2  to  10  nn  of  contaminant  buildup  can  cause  a 
degradation  by  a  factor  of  10  to  100  in  this  mirror's  scatter  performance 

Reference  25  reports  measurements  at  A  =  10.6  pm  of  BRDF  mirror 
scatter  increase  produced  by  ~  1  to  20-pm  cryodeposits  of  various  species 
Oxygen  deposits  consistently  caused  an  increase  in  scatter  with  increased 
deposit  thickness.  For  example,  the  data  indicate  that  an  O2  deposit  of 
1.7  pm  would  degrade  the  BRDF  to  the  CIRRIS  specification  of  1.8  x  10“^  sr_1 
at  2.3°.  In  contrast,  Nj  deposits  as  thick  as  ~  25  pm  did  not  ordinarily 
produce  any  measureable  BRDF  change.  In  a  few  cases,  however,  BRDF  increased 
by  ~  10  to  100.  This  was  presumably  associated  with  the  formation  of  an  N? 
frost,  as  opposed  to  that  of  a  transparent,  continuous,  homogeneous  smooth 
film. 

In  one  case,  the  scatter  from  an  N2  film  increased  markedly  after 
partial  evaporation,  perhaps  because  of  some  recrystallization  that  greatly 
increased  scatter  centers.  Subsequent  N2  deposition  onto  this  film  increased 
the  scattering  markedly.  This  worst-case  N2  behavior  would  indicate  that  the 
N2  cryodeposit  must  not  exceed  ~  0.05  un  (500  A)  in  order  to  meet  the  CIRRIS 
BRDF  requirement.  Because  better  information  is  not  available,  this  is  the 
safe  limit  recommended  here  for  contaminants  other  than  H2O  that  do  not  absorb 
appreciably.  It  is  probably  a  conservative  limit  because  it  is  based  upon  the 
worst  N2  behavior  observed.  Moreover,  it  is  also  subject  to  considerable 
uncertainty  because  it  is  arrived  at  frcm  an  assumption  that  BRDF  increase  is 
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proportional  to  cryodeposit  thickness  t  from  molecular  dimensions  to  dimen¬ 
sions  comparable  to  the  wavelength  X  of  the  incident  radiation.  In  reality, 
there  are  neither  data  nor  theories  to  support  such  a  linear,  or  other  par¬ 
ticular  functional,  extrapolation. 

Reference  25  also  reports  five  measurements  of  water  vapor 
depositions,  with  t  *  3  to  5  pm,  on  mirrors  at  temperatures  from  27  to  74  K. 
Three  of  these  measurements  revealed  no  increased  mirror  scatter  However, 
two  other  cases  revealed  a  mirror  scatter  increase  >  10^,  i.e.,  one  much 
larger  than  that  due  to  the  N2  frost.  A  linear  thickness  extrapolation  of 
this  behavior  gives  an  H2O  thickness  limit  of  ~  6  x  10-^  pm  (0.6  A)  before  the 
C1RRIS  BRDF  requirement  is  exceeded.  However,  this  linearly  extrapolated 
thickness  limit  seems  unrealistically  small,  which  indicates  the  failure  of  a 
linear  extrapolation  in  at  least  this  case.  Before  any  contamination,  the 
best  available,  state-of-the-art,  super-polished  mirrors,  with  BRDF  at  X  * 

10  pm,  considerably  better  than  the  CIRRIS  requirement,  have  root-mean- 
square  (RMS)  roughnesses  of  10  to  50  A.  Electromagnetic  boundary-value 
theoretical  expressions  given  in  Ref.  26  for  mirror  BRDF  can  be  used  to  esti¬ 
mate  RMS  surface  roughness  limits  for  which  should  be  highly  reflective 

at  wavelengths  where  it  absorbs  strongly.  These  theoretical  expressions 
depend  a  great  deal  on  the  surface  roughness  correlation  length  as  well.  They 
would  indicate  a  ~0.01  ym  surface  roughness  for  a  correlation  length  assumed 
to  be  much  less  than  a  wavelength  “0.1  ym  at  X  *  2.7  ym  (a  wavelength  near 
the  short -wave length  edge  of  the  CIRRIS  band  where  water  is  highly  absorbing) , 
in  order  to  meet  the  CIRRIS  BRDF  requirement  at  one  degree.  In  the  absence  of 
any  better  Information,  a  critical  thickness  limit  estimate  of  0.01  ym 
(100  A)  is  adopted  here  for  H2O  deposition  on  the  CIRRIS  primary  mirror. 


IV.  MONTE  CARLO  HELIUM-PURGE  EFFICIENCY  RESULTS 


A  summary  of  all  of  the  Monte  Carlo  calculations  is  given  in 
Table  2.  The  He  purge  efficiency  depends  on  a  number  of  variables,  viz  sensor 
geometry,  molecular  identity,  impact  speed  vc,  impact  direction  y,  He  mass 
flow  rate  m,  number  density  of  He  within  the  sensor  ,  and  temperatures  of 
the  He  effusion  sources  T^g  and  T  .  In  general,  these  results  show  the 
following  general  trends  in  transmission  of  contaminants  through  the  He  purge- 

(a)  t  decreases  monotonically  with  increasing  m  or  y  or  %e; 

(b)  x  increases  monotonically  with  increasing  vc  or  THg: 

(c)  x  (CIRRIS  IA)  >  t  (CIRRIS  I);  and 

(d)  x  (CO^)  >  x  (N^)  >  t  (H2O) ,  other  variables  being  equal. 

These  trends  are  discussed  in  more  detail  in  the  remainder  of  this  section. 

One  striking  feature  of  the  results  listed  in  Table  2  is  that  the 
purge  flow  of  He  provides  very  significant  protection  against  impinging  con¬ 
taminants  in  some  cases.  This  is  most  striking  for  CIRRIS  I.  In  many  of 
these  cases,  no  contaminants  reached  the  mirror  and  only  an  upper  limit  could 
be  set  on  t,  An  example  is  the  case  of  H2O  impinging  at  vc  =  6  x  10^  cm/sec, 
with  m  =  0.02  gHe/sec  and  T  <  0.002  (run  27).  Indeed,  1  <  0.001±0.001  for 
this  case  because  this  was  the  value  calculated  for  a  much  higher  impingement 
speed,  v  =  8  x  10"*  cm/sec  (run  26).  In  fact,  an  examination  of  wall  deposi¬ 
tion  numbers  in  the  different  calculations  indicates  that  x  <<<  0.001  for 
thermal  H20  at  m  =  0.02  g^e/sec,  perhaps  by  orders  of  magnitude.  Figure  6, 
for  example,  shows  results  for  H20  impinging  at  orbital  velocity.  Out  of  2000 
trajectories  followed,  one  reached  the  mirror,  whereas  371  others  deposited 
along  the  sensor  wall  with  an  approximate  exponential  decay.  In  contrast,  not 
a  single  one  out  of  a  total  of  1000  trajectories  even  entered  the  aperture  for 
thermal  H20  in  run  27.  Thus,  in  some  cases  the  He  purge  flow  does  indeed 
sweep  aside  completely  the  impinging  thermal  contaminants.  The  assumption 
discussed  earlier  of  a  unit  sticking  coefficient  is  clearly  inconsequential 
for  these  cases. 


The  remainder  of  the  discussion  presented  here  concentrates  on 
CIRRIS  IA,  the  instrument  that  will  be  used  on  an  upcoming  flight.  Both 
CIRRIS  I  and  CIRRIS  IA  were  modeled  as  telescopes  with  a  separation  T  =  94  cm 
between  the  sensor  aperture  and  the  primary  mirror.  Figure  7  compares  the 

CIRRIS  I  and  IA  entrance  apertures.  The  CIRRIS  IA  aperture  is  approximately 

2 

elliptical,  with  an  area  A  =  345  cm  ,  whereas  the  CIRRIS  I  aperture  was  cir- 

2 

cular,  with  A  =  222  cm  .  Because  the  Monte  Carlo  trajectory  program  required 
a  circular  aperture,  the  CIRRIS  IA  aperture  was  modeled  in  these  calculations 
as  circular,  with  a  diameter  equal  to  the  major  diameter  of  the  ellipse  and  an 
area  A  *  829  cm  .  At  a  given  flow  rate  m,  the  He  attenuates  incoming  contami¬ 
nants  more  efficiently  for  CIRRIS  I  than  for  CIRRIS  IA,  because  there  is  less 
sensor  aperture  area  to  protect  and  because  deflections  through  smaller 
angles,  on  average,  will  cause  the  contaminant  to  miss  the  sensor  aperture  or 
strike  the  inside  sensor  wall.  Thus,  the  CIRRIS  IA  results  presented  here, 
based  upon  the  much  larger  circular  aperture  area,  are  conservative  in  under¬ 
estimating  the  actual  attenuation  that  is  due  to  the  He  gas.  In  view  of  the 
values  of  the  three  areas,  in  fact,  the  results  presented  here  for  CIRRIS  IA 
should  be  very  conservative,  with  the  real  CIRRIS  IA  behavior  closer  to  that 
of  CIRRIS  I  than  are  the  computed  results  for  CIRRIS  IA. 

Figure  8  illustrates  some  of  these  remarks.  The  transmission  is  a 
strong  function  of  He  flow  rate.  The  He  provides  considerably  more  protection 
for  the  CIRRIS  I  geometry  than  for  the  CIRRIS  IA  geometry.  Even  for  CIRRIS 
IA,  however,  at  the  planned  flow  rates  of  0.02  g^g/sec,  the  He  flow  still  pro¬ 
vides  significant  protection  against  H2O  impinging  at  low  speed. 

Increasing  the  He  ejection  temperature  THg  should  increase  the  He 
momentum  while  decreasing  the  steady-state  density  in  front  of  the  sensor. 

The  effect  on  impinging  contaminants  can  be  seen  in  Table  2  by  comparing  runs 
27  versus  34,  32  versus  33,  and  28  versus  39.  Apparently  the  density  change 
Is  more  important.  In  all  cases,  THe  ■  20  K  was  considerably  more  effective 
than  TH@  *  273  K,  and  had  a  particularly  big  effect  on  the  N2  acceptance. 

Figure  9  Illustrates  that  H2O  transmission  shows  appreciable 
dependence  on  the  steady-state  density  of  He  within  the  sensor.  The  model  of 


Normal  View  of  Sensor  Apertures.  The  CIRRIS  I  aperture  is 
circular,  R  *  8.4  cm.  However,  the  CIRRIS  IA  elliptical 
aperture  was  modeled  by  a  larger  circular  aperture, 

R  *  16.25  cm. 


Transmission  of  H2O  Contaminant  as  a  Function  of  Steady-State 
He  Density  in  the  Sensor  for  Runs  44  through  47.  The  arrow 
indicates  that  NHe  \°  .55  x  10*  cm”  for  true  effusion  of  He 
out  of  the  ejector  ring. 


effusion  of  He  from  the  ejection  ring  into  the  sensor  actually  predicts 
N^e  ”0.55  x  10  cm  ,  so  that  the  true  value  is  likely  to  be  somewhere 
between  this  number  and  zero.  Contaminant  buildup  rates  that  are  estimated 
later  represent  the  conservative  analysis  where  N^e  is  taken  as  zero.  In 
reality,  the  buildup  of  ^0  should  be  significantly  slower  than  the  estimates 
for  H2O  provided  later  in  this  report,  as  a  result  of  a  finite  N^g.  The 
buildup  rate  of  N2 ,  however,  should  show  very  little  or  no  dependence  on  N^e, 
because  the  He  must  deflect  the  N2  before  it  enters  the  sensor.  There  might 
be  a  small  dependence  associated  with  the  small  increase  in  He  density  in 
front  of  the  sensor,  as  a  result  of  effusion  from  the  sensor  for  a  finite  N^e- 

The  transmission  of  H2O  is  plotted  versus  He  flow  rate  for  two 
speeds  of  contaminant  approach  in  Fig.  10.  This  illustrates  that  t  exhibits 
an  initial,  Beer's-law  exponential  decrease  with  increasing  m,  but  that 
significant  deviations  from  this  simple  dependence  appear  in  the  vc  =  6  x  10^ 
cm/sec  results,  where  the  attenuations  are  clearly  well  beyond  the  single¬ 
collision  regime.  Figure  11  illustrates  the  very  strong  dependence  of  t  on 
speed  of  approach  for  HjO.  This  is  reproduced  in  Fig.  12,  where  results  for 
N2  and  H2O  are  compared.  This  illustrates  that  the  contaminant  transmission 
is  more  sensitive  to  impact  speed  vc  than  to  molecular  identity.  For  this 
reason,  the  investigation  of  sensitivity  to  collision  cross  sections  was  not 
pursued  beyond  that  represented  by  runs  7,  48,  and  54.  Comparisons  of  runs  54 
and  33  are  interesting,  however.  Run  54  was  done  for  a  total  cross  section 
Independent  of  v,  with  a  value  equal  to  that  for  ^0  at  v  ”  8  x  10^  cm/sec  but 
with  an  Isotropic,  hard-sphere  differential  cross  section.  In  contrast,  run 
33  refers  to  the  correct  H2O  cross  sections,  so  that  the  total  cross  section 
increases  after  the  first  collision  has  slowed  the  ^0  down.  Nevertheless, 
run  54  produced  much  larger  attenuations  than  did  run  33  (x  *  0.12  and  0.43, 
respectively),  which  indicates  that  an  assumed  hard-sphere,  isotropic  differ¬ 
ential  cross  section  severely  overestimates  the  role  of  hard,  wide-angle 
collisions . 

As  noted  earlier,  ^0  is  expected  to  return  with  an  approximate 
thermal-speed  distribution  characteristic  of  the  outgas  source  for  8  >  90°. 

In  orbit,  Orbiter  surface  temperatures  can  vary  from  ~  200  to  375  K,  depending 
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Fig.  12.  Transmission  versus  v  for  HnO  and  N2 •  The  f^O  results  are 
the  same  as  In  Fig.  II;  the  Nj  results  are  from  runs  39,  49 
52,  and  53. 


on  solar  Illumination.  Figure  11  indicates  small  transmissions,  <0.01,  for 
average  thermal  return  speeds  corresponding  to  these  surface  temperatures.  In 
order  to  test  for  sensitivity  to  the  high-velocity  tail  of  the  Boltzman  speed 
distribution,  however,  the  T  (vc)  behavior  shown  in  Fig.  11  was  averaged  over 
a  thermal  speed  distribution  to  obtain  a  thermally  averaged  transmission  t. 

The  results  shown  in  Fig.  13  indicate  t  <  0.01  for  T  <  375  K. 

In  order  to  evaluate  Eqs.  (16)  and  (19),  t  and  n  were  evaluated  as 

a  function  of  y.  The  results  are  shown  in  Figs.  14  and  15  for  H20  and  N2 . 

Integration  of  the  t(y)  for  H00  at  v  =  8  x  10^  cm/sec,  shown  in  Fig.  14, 

-  c 

yielded 

bh2o-°-18  <21> 

The  calculations  of  r(y)  and  R  were  not  carried  out  at  lower  H20  impact  speeds 
because  of  the  very  large  corapution  time  that  would  have  been  required  to 
achieve  good  statistics.  Instead,  the  value  of  B  given  in  Eq.  (21)  was  used 
in  evaluating  Eqs.  (15a)  and  (17)  for  the  buildup  rate  of  H20  contaminant  film 
at  all  vc  values.  Again,  this  should  be  a  conservative  estimate,  since  B 
would  be  expected  to  decrease  with  decreasing  v  .  Runs  33  versus  67  and  34 
versus  68  in  Table  2  support  this  expectation,  although  the  statistics  are 
poor.  It  is  also  of  interest  to  compare  Eq .  (21)  with  a  rough  estimate  of 
2 

B  ■  A/L  .  For  the  CIRRIS  IA  model  employed  here,  this  rough  estimate  gives 
B  *  0.09,  whereas  it  gives  B  *  0.04  if  the  real  CIRRIS  IA  area  is  used. 

Again,  these  comparisons  serve  to  emphasize  that  the  results  calculated  here 
overestimate  conservatively  the  actual  contamination  buildup. 

The  N2  acceptance  results  shown  in  Fig.  15  indicate  a  much  weaker 
dependence  on  y.  Indeed,  the  He  purge  has  very  little  effect  on  acceptance  at 
the  highest  velocity  studied,  5.8  *  10^  cm/sec;  in  the  absence  of  He,  the 
behavior  would  have  been  1.0  cosy.  Acceptance  results  for  H2o  are  included  in 
Fig.  14  to  illustrate  that  this  behavior  in  Fig.  15  is  not  peculiar  to  N2 , 
Indeed,  the  results  in  Fig.  14  illustrate  that  He  does  little  to  prevent  fast- 
impinging  contaminants  from  entering  the  sensor,  but  that  it  does  serve  to 
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deflect  most  of  them  into  the  walls.  The  results  in  Fig.  15  clearly  show  that 
B'  of  Eq.  (19)  is  2 tt  at  vc  =  5.8  x  10^  cm/sec.  The  dashed  line  fit  to  the 
results  for  vc=  4.8  x  10^  cm/ sec  gives  B'  =  0.9.  A  functional  form  was 
assumed  to  join  these  two  limits’ 


B  =0.9  +  (5.4  cos  8) 
2 


where  8  is  calculated  from  Eq.  (12). 


V.  CONTAMINANT  BUILDUP  RATE  AND  CONCLUSIONS 


Figure  16  shows  the  buildup  of  H2O  computed  from  Eq.  (17)  with  the 
Orbiter  worst-case  1^0  environment  specified  in  Eq  (13).  For  a  total  H2O 
thickness  goal  of  <  ~0.01  pm  and  a  ~  100-hr  period  during  which  CIRRIS  IA 
will  be  uncapped  on  orbit,  growth  rates  of  t  <  10”^  ^/hr  should  be 
acceptable.  Thus,  Fig.  16  indicates  there  is  no  problem  associated  with  H20 
buildup  for  the  current  CIRRIS  IA  mission  plan:  m  =  0.02  gHe/sec,  8  =  90  to 
180°,  and  there  are  no  VCS  motor  firings.  Figure  16  also  illustrates  that  H20 
contamination  would  increase  significantly  if  the  sensor  were  pointing  into 
the  hemisphere  containing  the  incoming  wind.  The  second  curve  is  included  in 
Fig.  16  to  illustrate  what  part  of  this  increase  in  t  with  decreasing  8  is 
due  to  the  vc  dependence  of  t ,  rather  than  to  the  B  dependence  of  f^  ^ .  The 
functional  dependence  of  is  essenC^ally  unknown  except  at  0  =  0  and 

90°.  Thus,  the  true  t  probably  lies  between  the  two  curves  shown  in  Fig.  16, 
although  it  must  approach  the  curve  labeled  t  as  B  goes  to  0°.  Clearly  the 
CIRRIS  IA  mission  plan  is  wise;  it  would  be  unwise  to  plan  to  collect  all  data 
at  any  angle  8  <  90°.  If,  however,  it  might  be  desirable  to  collect  some 
limited  data  at  B  <  Fig.  16  can  be  used  to  estimate  the  corresponding 

contamination  penalty.  It  is  evident  that  some  data  collection  at  angles 
B  <  90°  would  be  acceptable. 

The  acceptance  n  of  N2  is  plotted  versus  vc  in  Fig.  17.  This  was 
used  in  conjunction  with  fN^(S)  fr°m  Eq .  (14)  and  ( B)  from  Eq .  (22)  to 
calculate  the  N2  steady-state  densities  (shown  in  Fig.  18)  within  CIRRIS  IA 
from  Eq.  (18).  For  a  CIRRIS  IA  primary  mirror  temperature  of  20  K,  condensa- 
tion  begins  at  ngg  *  5  *  10°  N2  molecules  cm”  ,  or  6  *  87°.  This  figure  also 
indicates  that  condensation  would  not  take  place  under  any  condition  at  a  wall 
temperature  greater  than  24  K.  Thus,  condensation  takes  place  only  at  the 
mirror  except  for  some  small  part  of  the  wall  near  the  mirror  which  may  have 
temperatures  of  20  to  24  K.  For  small  8  values,  the  N2  film  thickness  growth 
rate,  calculated  from  Eq .  (20),  is  shown  in  Fig.  19.  The  estimated-safe  N2 
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Fig.  17.  The  Acceptance  n  for  N2  versus  vc  for  Runs  39,  49,  52,  and  53 
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60  30° 

LINE-OF-S IGHT  ANGLE  FROM  WIND  S 


Steadv-State  N2  Density  within  CIRRIS  IA  versus  8  from  Eq .  (12) 
for  fN2(8)  from  E<1*  (TA)  and  B^2 ( 6)  from  Eq .  (22).  The  right 
vertical  axis  shows  the  maximum  temperature  at  which  the 
corresponding  number  density  shown  on  the  left  axis  would 
produce  condensation. 


0.8  1  2 

APPROACH  SPEED  vc.  105  cm/  sec 


the  CIRRIS  IA  Primary  Hirro 


thickness  discussed  earlier  was  0.05  pm,  so  that  t  should  be  less  than  5  x 
10-4  mn/hr  for  a  100-hr  mission.  Again,  the  mission  constraint  of  B  >  90°  is 
well-advised.  Figure  19  may  be  used  to  estimate  the  contamination  conse¬ 
quences  if  any  data  are  collected  at  smaller  B  angles. 

Buildup  rates  of  1^0  and  ^  contaminant  films  that  would  be 
produced  by  continuous  firing  of  a  VCS  engine  are  given  in  Table  3  for  the 
worst-case  return-flux  estimates  based  on  observations  reported  in  Ref.  17. 

As  discussed  earlier,  there  is  considerable  uncertainty,  given  a  large  range 
of  possible  return  speeds  and  f^O  and  N2  fluxes.  As  shown  in  Table  3,  this 
produces  a  large  range  of  uncertainty  in  possible  contaminant-film  buildup 
rates.  The  VCS  engines  are  used  when  the  Shuttle  Orbiter  is  not  in  the  stable 
gravity  gradient  mode;  they  are  operated  in  a  pulsed  mode  with  a  duty  factor 
much  less  than  unity.  The  corresponding  film  buildup  rate  is  obtained  by 
multiplying  the  entry  in  Table  3  by  D^.  For  a  typical  =  0.01,  the  thick¬ 
ness  buildup  rates  t  given  in  Table  3  may  be  read  as  total  buildup3  during  the 
100  hr  that  CIRRIS  IA  is  uncapped.  Even  these  worst-case  estimates  indicate 
that  H2O  buildup  is  not  a  problem,  but  that  N2  buildup  may  be  a  problem. 

A  more  careful  analysis  of  N2  buildup  requires  a  detailed  con¬ 
sideration  of  geometry,  because  vc  and  f  depend  on  the  location  and  thrust 
direction  of  the  VCS  engine  in  relation  to  the  Shuttle  Orbiter's  attitude  and 
the  sensor's  line  of  sight.  Table  4  provides  such  an  analysis  versus  B  for  a 
typical  case,  a  CIRRIS  IA  line  of  sight  straight  up  out  of  the  cargo  bay. 
Return  speeds  were  calculated  from  conservation  of  linear  momentum  and  energy 
for  an  N2  plume  molecule  that  returned  as  a  result  of  a  single  collision  with 
an  0  atom.  These  are  only  estimates.  In  reality,  a  spread  of  return  speeds 
would  be  produced  because  of  multiple  collision  effects  and  a  spread  in  the 
intersection  angle  of  the  collision  partners.  Combining  the  v£  values  with 
the  transient  N2  signals  reported  in  Ref.  17  gives  the  worst-case  fN^  and  t 
entries  in  Table  4.  Even  these  worst-case  entries  indicate  that  VCS  engines 


*It  should  be  noted  that  Nj  evaporation  while  the  motor  is  off  is  too  slow, 
given  the  20  K  temperature  of  the  primary  mirror,  to  be  of  any  consequence. 


Table  3.  Worst-Case  Estimates  of  H20  and  N2  Growth 
Rates  on  the  CIRRIS  IA  Primary  Mirror 
Due  to  the  25-lb-Thrust  VCS  Engine3 


Species 

VC> 

105  cm/ secb  1014 

f.  b 
cm”^  sr  *  sec-* 

Thickness  Buildup  t, 

pm/ hr 

h2o 

0.6 

0.4 

6  x  10-5 

h2o 

5 

4 

3  x  10-2 

n2 

0.5 

0.4 

2  x  10“3 

N2 

4 

3 

1 

^Thickness 

buildup  rate  refers  to 

continuous  engine 

burn : 

m  =  0.02  gHe/sec;  THg  =  273  K;  NRe  -  0 

^As  discussed  in  the  text,  this  is  the  estimated  range  of  possible 
return  speeds  and  worst-case  fluxes. 


Table  4  .  Best  Estimates  of  N2  Growth  Rates  on  the 
CIRRIS  IA  Primary  Mirror  Due  to  the 
25-lb-Thrust  VCS  Engine3 


f,  Molecules  Nj 
10^  cm-^  sr~*  sec-^- 


Thickness  Buildup  t, 
ixn/hr 


v  ,  Worst  Best 

1CH  cm/sec°  Casec  Estimate^ 


Worst 

Casee 


Estimate 


can  be  safely  fired  at  Df  =  0.01  during  the  entire  CIRRIS  IA  mission  so  long 
as  8  >  150°.  Moreover,  the  best-estimate  entries  in  Table  A  indicate  that 
contamination  buildup  due  to  VCS  engine  firings  at  a  typical  =  0.01  will  be 
inconsequential  so  long  as  the  sensor  is  looking  into  the  aft  hemisphere, 
i.e.,  if  0  >  90s.  These  best  estimates  were  calculated  by  extrapolating  the 
6  ■  0  to  80°  calculations  reported  in  Ref.  20  with  the  approximate  0  =  0  to 
180°  return  flux  angular  distribution  reported  in  Ref.  27.  Work  reported  in 
Ref.  28  indicates  that  even  this  should  overestimate  the  N2  return  flux  at 
these  large  0  angles,  because  N2  is  heavier  than  the  0  atom  collision  partner. 
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APPFNDIX 


Differential  cross  sections  have  been  calculated  by  the  method 
described  in  Ref.  3  for  the  scattering  of  He  atoms  by  H2O,  N?,  CO,  and  C02 • 
Results  are  presented  here  as  plots  of  o(9,v),  in  units  of  A  /sr,  as  a  func¬ 
tion  of  polar  scattering  angle  0.  Some  plots  are  presented  versus  collision 
energy,  E  =  pv  / 2 ,  rather  than  collision  speed.  For  completeness,  results  on 
He-0  are  also  included  here.  The  He-0  potential  function  given  here  is 
correct;  there  is  an  error  in  the  potential  function  that  appears  in  Ref.  3. 

Potential  Energy  Functions 


For  He  +  0(3P) : 


V(r)  *  tanh(ar-^)  VF(r)  +  tanh(br^)  VA(r) 


where : 


vf(r)  -  A«-6r 


vA(r)  =  e[-  f'  ,  exp  j-a[ (r/r  )-l] 
ct  ■“  d  in 


and  where 


a 

m 

102.0  bohr6 

b 

m 

0.001  bohr 

A 

X 

378  eV 

6 

m 

3.744  A-1 

e 
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1 
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*—< 

X 

00 

Csl 
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13.772 
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m 

3.08  A 
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For  He  +  CO  and  He  +  N-, : 

Both  potentials^  are  a  parametrized  form  called  the  Simons-Parr- 
Finlan-Dunham  (SPFD)  potential  given  by 


V(r) 


eLb0 


r  -  r  , 

tr-r-±y 


3 

I 

i=l 


r  -  r 


- “I1  -  1]:  *- 

r  J  r 


<  1.6 


The  parameters  are 


He  +  CO 
3.70  A 

2.28  x  10-3  eV 
39.2 

-  5.10 
11.31155006 

-  9.73209667 
1.52 


—  >  1.6 

r 

m 

He  +  N2 

r_  =*  3 . 65  A 

m 

e  =  2.27  x  10~3  eV 

bQ  =  31.9 

bl  =  -  4.72 
b2  -  10.05959840 
b3  -  -  8.52603029 
C6  =  2.10 


(b2  and  b3  are  determined  by  smoothness  conditions.) 


For  He  +  H2 : 

This  is  also  an  SPFD  potential7  with  the  following  parameters: 


r  -  3.38  A 

r 

e  -  2.27  x  10“3  eV 

bQ  -  31.4 

bj_  -  "  6.15 
b2  - 
bo  - 


17.6202007 

-18.5110554 


For  He  +  CO2 : 


This  potential  is  a  Morse-Spline-Van  der  Walls  (MSV)  potential' 
with  parameters  that  are  functions  of  the  angle  y,  as  follows: 

V(r,y)  =  c  exp[-6(r  -  rm)/rm]  {exp[-8(r  -  rm)/rm]  -  2}  (r  <  r}) 

v(r ,  y)  =  (r2  -  r)  [si(r2  -  r)2  +  S3  ]  +  (r  -  r^  [S2 (r  -  r^2  +  sj  (r2  <  r  < 

(r  >  r2) 


v,  .  C6  C8 
V(r,Y)  =  -  _  -  _ 

r  r 


where 

e  =  eo  +  e2  p2^C0SY) 

rm  =  rm0  +  rm2P2^C0SY) 

C6  =  C6(0)  +  C6(2)P2(cosy) 

C8  =  c8(°)  +  C8(2)P2(cosy) 

The  function  P2(cosy)  is  the  second-order  Legendre  polynomial  and  the  values 
of  the  parameters  are 

Eq  =  1.05  meV 

e2  =  -  0.67  meV 


rm0 

3. 

.84 

A 

rm2  = 

0, 

.99 

A 

C6(0)  - 

9, 

.98 

eV 

A6 

C6(2)  * 

2. 

.31 

eV 

A6 

C8(0)  - 

46. 

.4 

eV 

A® 

Cg(2)  - 

48. 

,4 

ev 

A® 

4.59 


(1  + 


ln(2) 


)  r„ 


1.6  r„ 


r2 ) 
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The  spline  parameters  ,  S2 ,  S3,  and  are  determined  by  smoothness  condi 
tions  at  r^  and  i.e.,  the  potential  and  its  first  derivative  are 

continuous  at  these  points. 
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LABORATO  OPERATIONS 


The  Aerospace  Corporation  functions  as  an  "architect -engineer"  for 
national  security  projects,  specializing  in  advanced  military  space  systems. 
Providing  research  support,  the  corporation’s  Laboratory  Operations  conducts 
experimental  and  theoretical  investigations  that  focus  on  the  application  of 
scientific  and  technical  advances  to  such  systems.  Vital  to  the  success  of 
these  investigations  is  the  technical  staff's  wide-ranging  expertise  and  its 
ability  to  stay  current  with  new  developments.  This  expertise  is  enhanced  by 
a  research  program  aimed  at  dealing  with  the  many  problems  associated  with 
rapidly  evolving  space  systems.  Contributing  their  capabilities  to  the 
research  effort  are  these  individual  laboratories: 

Aerophysics  Laboratory:  Launch  vehicle  and  reentry  fluid  mechanics,  heat 
transfer  and  flight  dynamics;  chemical  and  electric  propulsion,  propellant 
chemistry,  chemical  dynamics,  environmental  chemistry,  trace  detection; 
spacecraft  structural  mechanics,  contamination,  thermal  and  structural 
control;  high  temperature  thermomechanics,  gas  kinetics  and  radiation;  cw  and 
pulsed  chemical  and  excimer  laser  development  including  chemical  kinetics, 
spectroscopy,  optical  resonators,  beam  control,  atmospheric  propagation,  laser 
effects  and  countermeasures. 

Chemistry  and  Physics  Laboratory:  Atmospheric  chemical  reactions, 
atmospheric  optics,  light  scattering,  state-specific  chemical  reactions  and 
radiative  signatures  of  missile  plumes,  sensor  out -of -f leld-of -view  rejection, 
applied  laser  spectroscopy,  laser  chemistry,  laser  optoelectronics,  solar  cell 
physics,  battery  electrochemistry,  space  vacuum  and  radiation  effects  on 
materials,  lubrication  and  surface  phenomena,  thermionic  emission,  photo¬ 
sensitive  materials  and  detectors,  atomic  frequency  standards,  and 
environmental  chemistry. 

Computer  Science  Laboratory:  Program  verification,  program  translation, 
performance-sensitive  system  design,  distributed  architectures  for  spaceborne 
computers,  fault-tolerant  computer  systems,  artificial  Intelligence,  micro¬ 
electronics  applications,  communication  protocols,  and  computer  security. 

Electronics  Research  Laboratory:  Mi croelectronics ,  solid-state  device 
physics,  compound  semiconductors,  radiation  hardening;  electro-optics,  quantum 
electronics,  solid-state  lasers,  optical  propagation  and  communicat ions ; 
microwave  semiconductor  devices,  microwave/millimeter  wave  measurements, 
diagnostics  and  radiometry,  microwave/millimeter  wave  thermionic  devices; 
atomic  time  and  frequency  standards;  antennas,  rf  systems,  electromagnetic 
propagation  phenomena,  space  communication  systems. 

Materials  Sciences  Laboratory:  Development  of  new  materials:  metals, 
alloys,  ceramics,  polymers  and  their  composites,  and  new  forms  of  carbon;  non¬ 
destructive  evaluation,  component  failure  analysis  and  reliability;  fracture 
mechanics  and  stress  corrosion;  analysis  and  evaluation  of  materials  at 
cryogenic  and  elevated  temperatures  as  well  as  in  space  and  enemy-induced 
environments. 

Space  Sciences  Laboratory:  Magnetospher ic ,  auroral  and  cosmic  ray 
physics,  wave-part ic le  interact  ions ,  magnetospherlc  plasma  waves;  atmospheric 
and  ionospheric  physics,  density  and  composition  of  the  upper  atmosphere, 
remote  sensing  using  atmospheric  radiation;  solar  physics,  infrared  astronomv. 
Infrared  signature  analysis;  effects  of  solar  activity,  magnetic  storms  and 
nuclear  explosions  on  the  earth's  atmosphere,  ionosphere  and  magnetosphere; 
effects  of  electromagnet ic  and  particulate  radiations  on  space  systems;  space 
instrumental  Ion. 


